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Following the consensus of a workshop in Turbulence Modeling for Liquid Rocket 
Thrust Chambers, the current effort was undertaken to study the effects of second-order 
closure on the predictions of thermochemical flow fields. To reduce the instability and 
computational intensity of the full second-order Reynolds Stress Model, an Algebraic 
Stress Model (ASM) coupled with a two-layer near wall treatment was developed. Various 
test problems, including the compressible boundary layer with adiabatic and cooled walls, 
recirculating flows, swirling flows and the entire SSME nozzle flow were studied to assess 
the performance of the current model. Detailed calculations for the SSME exit wall flow 
around the nozzle manifold were executed. As to the overall flow predictions, the ASM 
removes another assumption for appropriate comparison with experimental data, to account 
for the non-isotropic turbulence effects. 
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— Efficiency and stability 



APPROACH 
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APPROACH -Pressure-Strain Term 
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Two-Layer Wall Treatment 




METHODOLOGIES 
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Time Marching 



Validations 
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Compressible Flat Plate Flow 
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Fig. 1 Variation of Cf/C/o with for adiabatic wall boundary condition. 




Compressible Flat Plate Flow 
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Fig. 2 Variation of Cf/Cfu with for M 











Backward Facing Step 
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Fig. 5 Reynolds stress profiles for the backward-facing step turbulent 
flow (9:1), with data from (31] 





Confined Swirling Flows 
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Fig. 6 Decay of mean axial centerline velocity. 



Confined Swirling Flows 
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Fig. 8a Radial profiles of turbulent intensity(\/u^) 



Confined Swirling Flows 
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Fig. 8b Radial profiles of Reynolds stressfnV) 
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of temperature for ASM with two-layer model. 
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8 - STEP REACTIONS 




A 

N 

E 

M + 0 2 

===== 0 + 0 

0.72000E + 19 

-1.0000 

117908 

m + h 2 

==== H + H 

0.55000E+19 

-1.0000 

103298 

m + h 2 o 

==== H +0H 

0.52000E+22 

-1.5000 

118000 

O + H 

==== OH 

0.71000E+19 

-1.0000 

0. 

H 2 0 + OH 

==== H 2 0 + H 

0.58000E+14 

0.0000 

18000 

H 2 + OH 

==== H 2 0 + H 

0.20000E+14 

0.0000 

5166 

o 2 + h 

==== OH + O 

0.22000E+15 

0.0000 

16800 

h 9 +o 

==== OH +H 

0.75000E+14 

0.0000 

11099 


k = ATN ex p(-% T ) 

with kin cm3- mole-1 . S -1 and Ein cal • mole-1 


From R.C. Rogers and Chinitz, ' Using a global hydrogen-air model 

in turbulent flow calcultions ', AIAA J. , vol. 21, pp. 586-592, 
1983 . 
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SSME Nozzles 
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Figure A. 5 Sample SSME Nozzle Flow Inputs and Results Turbulent, 

8-Step Kinetics. 

ISP - 452.78 sec. Exp. ISP=453.3sec 


1738 






nozzle exit. 
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distance from sub-domain inlet along the wall, inches 


Fig. 17(a), Pressure levels along the wall near the nozzle exit using the 
ASM and k - e models 
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Fig. 17(b), Effects of wall temperature on the wall pressure 
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Summaries 
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3-D Calculations Are Desirable 
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